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Abstract. We analyze atom-surface magnetic interactions on atom chips where the magnetic trapping 
potentials are produced by current carrying wires made of electrically anisotropic materials. We discuss a 
theory for time dependent fluctuations of the magnetic potential, arising from thermal noise originating 
from the surface. It is shown that using materials with a large electrical anisotropy results in a considerable 
reduction of heating and decoherence rates of ultra-cold atoms trapped near the surface, of up to several 
orders of magnitude. The trap loss rate due to spin flips is expected to be significantly reduced upon 
cooling the surface to low temperatures. In addition, the electrical anisotropy significantly suppresses the 
amplitude of static spatial potential corrugations due to current scattering within imperfect wires. Also the 
shape of the corrugation pattern depends on the electrical anisotropy: the preferred angle of the scattered 
current wave fronts can be varied over a wide range. Materials, fabrication, and experimental issues are 
discussed, and specific candidate materials are suggested. 



PACS. 37.10.Gh Atom traps and guides - 39.25.-|-k Atom mainpulation 
in metals and alloys - 03.75.-b Matter waves 



72.15.-v Electronic conduction 



1 Introduction 

The use of atom chips for trapping, cooling, manipulat- 
ing and measuring ultra-cold atoms near surfaces has at- 
tracted much attention in recent years [HI]. The mono- 
lithically integrated micro-structures on the chip lead to 
tight potentials, which can be tailored over length scales 
comparable to the atoms' de-Broglie wavelength. This has 
simplified the production of ultra-cold quantum degener- 
ate gases [3l|4] , enabled high precision experiments such as 
atom interferometry [SUS], and has established ultra-cold 
atoms as a probe of the nearby surface of the atom chip [7] . 
To date, atom chip traps are mostly magneto-static. How- 
ever the use of electrostatic, radiofrequency (rf), or light 
potentials for atom manipulation on atom chips is also 
evolving rapidly (e.g. [SlIWTOlinifnilT^ ) . The chip plat- 
form is also considered a candidate for the development of 
quantum technologies in the field of quantum information 
processing and communication, as well as in interferome- 
try based sensors. The advantages of being in close prox- 
imity to the surface are hindered, however, by harmful 
processes originating from the surface itself, which cur- 
rently limit the capabilities of this platform [l4jfT5] . 

Random motion of thermal electrons within the nearby 
surface (known as Johnson or thermal noise) leads to mag- 
netic field fluctuations, affecting the trapping potential 
and resulting in trap loss, heating, and decoherence. As 
a general rule, the thermal noise increases as the atoms 
approach the surface; therefore it sets limits to the con- 



finement of micro-traps and to the coherence time in sen- 
sitive interference experiments [Tl[2] . The issue of trap loss 
due to spin flips has received the largest attention so far. 
Several studies were done on characterizing the lifetime 
of trapped atoms at different atom-surface distances, and 
with different material surfaces [T6l[T7l[T8plp0pT|[22l[23] . 
and the theory was found to fit well to experimental data. 
Recently, some theoretical work studied the reduction of 
trap loss, when the surface is either made of certain metal 
alloys [22 or superconductors [24 l [25 l [26 l [27 p 8 l [29] . while 
cooling the surface to cryogenic temperatures. Very few 
experiments were done on surface noise induced decoher- 
ence 



Current scattering (change of current flow direction) 
from surface or edge roughness or from bulk inhomo- 
gcneitics leads to static spatial corrugations of the mag- 
netic trapping potential [151131] . The consequent atomic 
density variation results in the fragmentation of the atom 
cloud into several separated components as it is brought 
closer to the surface. Much work has been done theoreti- 
cally and experimentally in order to understand the origin 
of fragmentation P^[?T1I5^1I331IM] . Improved fabrication 
has led to a considerable reduction of the corrugations, 
and schemes have been suggested for further improvement 
|151l35j . Recently, highly organized corrugation patterns 
were observed and analyzed [34ll36| . when a ID cloud was 
scanned across lithography patterned wires of different 
characteristics. A prominent feature, observed for different 
types of wires, were long range wave fronts along which 
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current scattering is maximal, oriented at ±45° to the wire 
axis. This was explained in terms of maximally scattering 
Fourier components of random conductivity fluctuations 
in the wire [341136]. 

In this paper we study both thermal noise and static 
potential corrugation in the context of atom chip conduc- 
tors that are electrically anisotropic. These have different 
conductance along the different crystal axes, contrary to 
the isotropic conductivity of normal metals. The prop- 
erties of such materials, including ruthenates, cuprates, 
graphites and others, have been studied thoroughly (see 
Table 1 below). Regarding thermal noise, we show here 
that electrically anisotropic materials lead to a significant 
reduction of the heating and decoherence rates. Trap life- 
times can be increased by cooling the chip surface to low 
temperatures. With respect to potential corrugation, we 
show that the orientation of current scattering wave fronts 
within wires can be tuned by the electrical anisotropy. For 
high anisotropy, fragmentation is expected to be strongly 
suppressed. 

In section [2] we discuss the implications of using 
anisotropic materials on the thermal noise produced by 
a surface, and analyze the resulting trap lifetime, heating, 
and decoherence processes. We then extend in section [3] 
the formulation used previously [341136] to analyze poten- 
tial corrugation and fragmentation to the case of electri- 
cally anisotropic materials. Finally, in section[4]we discuss 
issues of materials and fabrication, important for the de- 
sign of experiments to test the theory, as well as for future 
engineering of quantum technology devices. 

2 Loss, heating, and decoherence due to 
thermal noise 

The theory describing the affect of thermal noise on ultra 
cold atoms is now well established for isotropic materials 
(e.g. [TlfT4 l l37 l l38 l l39] ) . Other models, based for example 
on non-local electrodynamics, have also been developed 
[40] , In this section we derive and discuss a generalized 
theory of thermal noise coupled to magnetically trapped 
atoms, which includes the case of electrically anisotropic 
materials. Our method incoherently adds up elementary 
current sources in the material, within the quasi-static 
approximation [37] . Although approximate (up to correc- 
tions on the order of a factor 2), this approach can be used 
to obtain closed- form solutions for many wire geometries. 
This is in contrast to methods based on Green functions, 
which are accurate but limited to simple geometries such 
as a half-space or a laterally infinite layer. It should be 
noted that the quasi-static approximation is valid when 
the skin depth S = y^2/(TofiouJ (where ao, no, and ui are 
the electrical conductivity, permeability of free space, and 
the radiation frequency, respectively) is much larger than 
both the atom-surface distance d and the wire thickness 
H. At room-temperature the skin depth of Au at 1 MHz 
is on the order of lOfim. Electrically anisotropic materi- 
als usually have a much higher resistivity, and their skin 
depths are orders of magnitude larger at the same fre- 
quencies (e.g. for graphite S w 1mm, see Table 1 below). 



The quasi-static approximation thus applies even better 
to the latter. This approximation has been corroborated 
to a high degree of accuracy in experiment [23t[2l1 and 
theory [iT] . 

2.1 Physical processes 

We now introduce the different physical processes result- 
ing from the coupling of thermal noise to the trapped 
atoms. We show that the important quantity is the noise 
power spectrum, and calculate the latter in the follow- 
ing section. The thermal radiation couples to the atoms' 
magnetic moment by the Zeeman interaction V {x,t) = 
—fi ■ B (x, t), where fi is the magnetic dipole moment op- 
erator. This operator can be written for convenience as 
fj, = fisgpF^ with fiB Bohr's magneton, gp the Landc 
factor of the appropriate hyperfine level, and F the to- 
tal spin operator. As the magnetic field noise is ran- 
dom in time and space it averages to zero, and its ef- 
fect is expressed through its two-point correlation func- 
tion {V {xi,t) V {x2,t')), between different points in time 
t,t' and space xi,X2. 

The field and magnetic moment being vectors, one ac- 
tually needs the cross-correlations between their compo- 
nents to characterize the different processes induced by 
the noise (spin flips, heating, and decoherence). It follows 
from time-dependent perturbation theory that the rate for 
a transition from an initial state |0) to a final state |/) of 
the system is given by [42] 

ro-.f = ^J2 /d(t-Oe'-°^'(*-*')x (1) 
hi 

X ((0|/i,B,(x,t)|/)(/|A*ji?,(x,t')|0)) 

where the indices i,j label the Cartesian components of 
the magnetic field. The argument x of the magnetic field 
in Eq. ([T]) is the atomic position operator. The matrix el- 
ements thus involve spatial (overlap) integrals over the 
spatial wave function and magnetic field (e.g. see Eq. 
([7]) below). These integrals can be written in terms of 
the two-point correlation function of the magnetic field. 
The integral of this correlation function over the time 
difference {t — t') is related to the cross-spectral density 
5^(xi,X2; Wo/) of the magnetic fluctuations at positions 
Xi , X2 , and at the transition frequency coq / [T41I43] . This 
is defined as 

S'j^{xuxr,cj)= (2) 

/oo 
d{t-t')e-"^^'-''^ {B,{x,,t)B,{x2,t')), 
-OO 

or equivalently by 

{B* {xi,iu)Bj {X2,LU')) = 27rS{uj^Lo') S'i {xi,X2]C0) . 

(3) 

In cases where the spatial degrees of freedom can be 
considered classically, x in Eq. ^ can be taken as the co- 
ordinate of the atoms, giving rise to a position-dependent 
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transition rate between spin states. Instead of the wave 
function overlap integrals, the coordinates are then taken 
as Xi = X2. For calculating at the trap center one takes 
Xi ^ X2 = r [44] . 

The noise power spectrum S'^(xi,X2; wq/) is quite flat 
at all relevant frequencies which correspond for example 
to transitions between Zeeman magnetic sub- levels (rf), 
or the trap vibrational states (Hz to kHz range). Hence 
a low- frequency limit wq/ ^ can be taken. If this spec- 
trum is not flat in the relevant frequency range, the time 
dependence of loss and decoherence processes is more com- 
plicated and does not reduce to a simple rate [H] . 

Spin flips are transitions from trapped (|0)) to un- 
trappcd (I/)) internal states. These states are eigenstates 
of the spin operator component parallel to the quantiza- 
tion axis, defined by the static trapping field at the bottom 
of the trap. Due to the form of the dipole moment opera- 
tor fii, the Zeeman interaction can induce a spin flip only 
when the direction of the magnetic field fluctuation is per- 
pendicular to the quantization axis. Hence, in this case, 
we can rewrite ([Ij as 

7.pinmp= E ^5^"(r;c.o/), (4) 

where we sum over the perpendicular components /, m 
only, and take the matrix elements fii = fiBgF{0\Fi\f) 
of the dipole moment between the states |0) and |/). The 
spin flip rate can be measured from the lifetime of a mag- 
netic micro-trap, and by varying .the trap distance, one 
can discriminate loss due to surface induced magnetic field 
fluctuations, against trap loss of different origin [T8 l [T9 l l20 l 

mm- 

Decoherence of a quantum superposition state can oc- 
cur without changing the occupation of the states involved 
in the superposition, affecting only the phase. We distin- 
guish between spin decoherence and spatial decoherence. 
The former involves the change of relative phase of two 
internal states and |2)) in a superposition at the same 
spatial location r, while the latter involves the change 
of relative phase of two spatially separated components 
of the atom cloud (positions Xi and X2), trapped in the 
same internal state. Here we consider only 'classical deco- 
herence' in which the phase change arises from the fluc- 
tuations in the Zeeman shift. Such shifts occur, to lowest 
order, for magnetic field fluctuations parallel to the atom's 
magnetic dipole moment. Hence for both types of decoher- 
ence processes we need only the parallel component of the 
same power spectrum appearing in the spin flip case (Eq. 
([l])). The spin decoherence rate can thus be written as 

'Yspin decoherence — '2H'^ ' 

with = (2 |/i|| I 2) — (1 I 1) being the differential 
magnetic moment of the two trapped states, and (r; 0) 
the low-frequency limit parallel component of the noise 
spectrum, at the trap center r. This decoherence rate can 
be measured from the reduced contrast of interference 



fringes in a series of Ramsey spectroscopy experiments, 
performed as a function of the atom-surface distance. 

For the case of a spatially separated superposition 
state, the rate of decoherence of the relative phase of the 
atomic states between two points Xi,X2 involves the cor- 
relation function of the difference in the magnetic fields 
B\^{xi,t) — B\^{x2,t'), and can be written as 

_ t^l 

Tspatial decoherence — [Sll + S22 ~ 26*12] , (6) 

where we denote for convenience Sij = S\\{xi,Xj;uj —^ 0), 
and assume a correlation spectrum symmetric in Xi, X2 
and flat in frequency. We again take only the field com- 
ponents parallel to the quantization axis, that shift the 
relative phase between the two parts of the wave func- 
tion, /^ii is the magnetic moment matrix element of the 
single internal state. The low-frequency limit is valid here 
as the "scattering cross section" for transitions involved 
in the decoherence process is in practice independent of 
the frequency p^[39] . 

This decoherence rate can be measured by studying the 
interference pattern contrast reduction in double-well ex- 
periments, where the atom cloud is separated in two parts, 
which arc then held at a fixed separation for a given time. 
The interference pattern is obtained by overlapping the 
cloud parts upon release from the trap. The interference 
contrast is directly related to the product of split time 
and decoherence rate; it can be measured as a function of 
the atom-surface distance by repeating the measurement 
at different heights. 

Finally, heating of the trapped atoms (as a result of 
exciting external degrees of freedom while retaining the 
same internal state) can also be caused by the coupling 
to the thermal radiation. The transition rate of atoms ini- 
tially in the ground vibrational state |0) to higher states 
/) with energy splitting Tiwof is of the form [14] 

fJ-l f 

^"^/ / dxidx2MjQ{xi)Mfo{x2)Si\ {xi,X2;ujof) , 

(7) 

where we find once more the spin operator matrix ele- 
ment in the direction parallel to the quantization axis, 
and now also the wave functions of the levels involved in 
the transition M/o(x) = '4'f{x)il;o{x). The spatial integra- 
tion here provides a probe of the spatial correlation of the 
magnetic field noise. In practice, it is enough to consider 
transitions from the ground state to cither of the first two 
excited levels [l]. 

Thus, we see that the important term common to all 
rates is the spectral density of the magnetic field fluctu- 
ations or power spectrum 5^ {xi,X2;uj). This quantity 
holds all of the relevant information about the magnetic 
field fluctuations leading to the harmful processes, while 
the other terms in each of the rates describe the coupling 
of the noise to the atoms through the magnetic dipole mo- 
ment. Furthermore, we see that a measurement of either 
the spin or spatial decoherence rates, or of the heating 
rate, will give strong indications to any of the other two 
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of these three processes, as they ah depend on the noise 
power of the same field component. 



2.2 Calculation of the noise power spectrum 

As presented in the previous section, the noise power spec- 
trum is related to the cross correlation function of the 
magnetic field fluctuations (Eqs. ([2|),([3])). Before calculat- 
ing this correlation function, we define the coordinate sys- 
tem such that the wire length L is along the x direction, 
its width W along the y direction, and its thickness H 
along z. As in typical magnetic traps a bias field polar- 
izes the atoms along the wire, the quantization axis in our 
analysis is along the x axis. 

We can make a distinction between two types of 
anisotropic materials. Materials with a 'layered conduc- 
tance' are relatively good conductors along two axes and 
have one axis of bad conductance: ayy <C a-xx ctzz- We 
always assume that the wire is aligned to one axis of good 
conductance, so that current may flow easily and create a 
magnetic trap. Materials that have only one direction of 
good conductance, ayy ~ a^z < cf^x or a^z < cTyy < o^xk, 
will be denoted as 'quasi-lD conductors'. Specific materi- 
als are discussed in Sect. 14.1] and Table [TJ 

Turning now to the magnetic field fluctuations, we use 
the expression for the current cross correlation function 
03116], 

AirheoLo^n [lo) 5 {ijJ — uo') Ime^ {xi^w) 5 [xi — X2) , (8) 

where n{Lj) — l/(e''B'^ — 1) is the Bose-Einstein occupa- 
tion number, and where the dielectric tensor 



eouj 



(9) 














(10) 










is proportional to the conductivity tensor for homogeneous 
media. When the anisotropic crystal axes are aligned with 
the wire, the latter is diagonal. 



This property will be used below. We assume here that 
the current in the material responds locally to the elec- 
tric field (Ohm's law), ji(x) ~ aij{x)Ej{x.). The vector 
potential and its correlation function in the quasi-static 
approximation are given by 

A{x,u;) = ^Jdx''^ (11) 

{A* {x,,^)A, {X2,u')) cc /^da^'^^^^Jg^. (12) 

Here Xi, Xi denote the two spatial locations for which we 
take the correlation function, whereas x' is the integra- 
tion variable, such that we sum the contribution from all 



points within the material volume Y , which are at dis- 
tances \xi — x'\ (i = 1, 2) from the locations a^i, X2- Note 
that this formula neglects jumps in A due to surface cur- 
rents at the metal-vacuum interface. In order to calculate 
the correlation function of the magnetic field fluctuations, 
we take the curl of the vector potential correlation func- 
tion, once with respect to X\ and once with respect to X2. 
Writing this in tensor form we get. 



(B* {x-Y,Uj)Bj (X2,J)) CX 
- / dx'eilm^jnpdl,ld2,n-, T, 

2 J ' ' ' \xi-x'\ 



\X2 - X' 



Hy(13) 



using the Levi-Civita symbol eijk and summing over re- 
peated indices. We defined the integral holding the con- 
ductivity tensor and the geometry terms as Bij for conve- 
nience, and the symbol da^i (a 1, 2) means a derivative 
with respect to x^ in the direction of its /th-component. 
Performing the derivatives we obtain 



B, 



(14) 



where the geometry of the system enters through the 
quantity Xin (correcting a missing factor 1/2 in Eq. (A. 6) 
of Ref. [37]): 



X,, 



dx 



,ixi-x')^ {X2-X')^ 



\xi - X'\^ \X2 



(15) 



We assume here a homogeneous medium (i.e., the compo- 
nents ijij are spatially constant within the material volume 
V). 

Again limiting the discussion to the 'aligned' case pU]) . 
Eq. (fT4)) is simplified, and in fact for every pair of i,j 
we need to sum only two integrals. Considering the wire 
geometry to be such that the atoms are located above the 
center of a very long wire (L ^ H,d), the only non-zero 
elements are Bu (i = x,y,z) due to symmetry. Hence we 
obtain 



Bxx — ^zzXy 



CTyyXzz 



Bzz '^yyXxx '^xxXyy 



and it is convenient to define 



^ij — Bij I CTxx-i 



(16) 



(17) 



assuming the good conductivity to be along x. 

We see that in contrast to the isotropic case, where 
one has a single conductivity ctq for all field components 
-By , here the conductivities an give different weights to 
the geometry-dependent factors Xa. 

We can now write the full expression for the power 
spectrum, using ^ and collecting all of the prefactors 
neglected in dEJ) , ([13]) as 



{Xi,X2 



-ImexxYij, (18) 

Anoj/c 
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where, following |37| , wc have normalized the power spec- 
trum to Planck's blackbody formula 



(19) 



As the relevant frequencies are low, the high temper- 
ature limit of the Planck function is applicable, and thus 
the expression for the power spectrum approaches 



5*5 {xi,X2;uj) 



47r2e2c4 



(20) 



This expression has the same form as the result for the 
isotropic case [37|, however here the tensor Yij holds the 
anisotropic terms Bij in it. 



2.3 Trap loss due to spin flips 

Examining Eq. wc sec that B^x can be considerably 
reduced if it involves two badly conducting axes, Uyy^Ozz 
being much smaller than Uxx- This does not lead, however, 
to a reduction of the spin flip rate, as the Bxx noise is 
parallel to the quantization axis. 

For the two perpendicular components, Byy and Bzz, 
we find that both mix the highly conducting Uxx and the 
low conductivity terms. The geometrical factors Xij have 
been analyzed in detail (see appendix A in [23]) for the 
case of rectangular wires. From this analysis it emerges 
that for any reasonable wire geometry, all of the non-zero 
Xij factors are on the same order. Thus the main differ- 
ence in the noise components is due to the difference in 
conductivity terms, where conductivity Uxx is dominant. 
Consequently, the improvement to the trap lifetime using 
anisotropic materials at room temperature is expected to 
be at most on the order of ~2. 

One could, of course, rotate the anisotropic crystal by 
90° such that the better conducting axis is perpendicu- 
lar to the wire. This may significantly reduce the spin flip 
rate for some materials, as now the conductivity along 
the quantization axis (along the wire) is smaller. However, 
this is not practical if the wire is to be used as a current 
carrying structure creating the magnetic trapping poten- 
tial, since the power consumption will be signiflcantly in- 
creased. This and other considerations on materials and 
experimental design will be discussed in section |4l 

The trap lifetime may be improved nevertheless by 
cooling the anisotropic material to cryogenic tempera- 
tures, as has been shown also for certain metal alloys [23] . 
The resulting improvement in lifetime is shown in Fig. [U 
where we compare as an example the lifetime due to ther- 
mal noise from an SrNb03.4i wire (a quasi-lD material, 
see Table 1 [55]) to that from a Au and an Ag:Au alloy 
wires. We see that for SrNb03.4i an improvement of two 
orders of magnitude in lifetime is expected upon cooling. 
This is due to the behavior of the product T ■ Gxx (T) (Eq. 
(|20)) ') for this and similar materials. Wc note that although 
at low temperatures the skin depth is decreased due to the 
smaller resistivity of the material, the quasi-static approx- 
imation is still valid. For a detailed discussion see [23]. 




100 ISO 200 

T [K] 

Fig. 1. Improved trap lifetime upon cooling of the surface. 
Comparison of a standard Au wire (dahsed blue) with wires of 
similar geometry made of an Ag:Au alloy (dashed-dotted red) 
[23] and the quasi-lD SrNb03.4i (solid black). The long life- 
time for the anisotropic wire at room temperature is due to its 
rather high residual (low-temperature) resistivity even along 
its best axis (oriented along the wire). See section IT2l for the 
implications of highly resistive wires in realistic experiments. 
No limiting lifetime level due to collisions with background 
gas atoms is included. Losses due to this process are expected 
to be less significant in cryogenic experiments as typically the 
background pressure is such experiments is smaller. Wire di- 
mensions are: width W = lO/xm, thickness H = 2.15/xm, and 
atom-surface distance d — 5/im. 



2.4 Heating and decoherence 

In contrast to the case of spin flips, the decoherence and 
heating rates, Eqs. ()5l6|7p . depend only on the noise in the 
parallel held component Bxx- As shown in section [^31 this 
component may be strongly reduced for highly anisotropic 
materials. To quantify this, consider the ratio 



paniso _ V 

^xx _ "zz^yy 

Qiso ~ 

XX 



(TyyXzz 



^XX {Xyy Xzz) 



(21) 



where the reference level is an isotropic conductivity set 
to the 'good axis', ctq = cJxx- The possible improvements 
depend on the relative magnitudes of the anisotropic 
conductivities. Four cases can be distinguished, as il- 
lustrated in Fig. [2] For materials with layered conduc- 
tance, the worst choice is to have the second good con- 
ducting axis in the chip plane, along the wire's width: 
'^vv ~ '^xx 3* cTzz (dashed red line). The ratio (|2ip 
then tends to (1 -f Xyy/Xzz) ^ which is not signiflcantly 
smaller than unity and where the conductivity anisotropy 
<^xx/<^yy actually does not enter. With the other choice, 
having the badly conducting axis along the wire width 
(dashed-dotted blue), we get a reduction of about one 
order of magnitude for a small wire geometry. Materi- 
als that are quasi-lD conductive have a much larger po- 
tential for noise suppression: for comparable 'bad axes', 
'^yy ^ '^zz ^ <^xx^ Eq- (|2ip scales inversely with the 
anisotropy ratio r = (Jxxl<^yy which may be very large. 
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lo" lo' 10^ ^ lo' 10^ lO' 



Fig. 2. Lines: Spin decoherence rate 7spin dccohorcnco as a func- 
tion of electrical anisotropy r = axx/cfyy, for layered and quasi- 
ID conducting materials. Wire parameters were d = 5^m, 
W = lO/xm, H = 2.15pim, and surface temperature of T = 
SOO-R'. For these lines, the good conductivity along the wire 
was assumed to be identical to that of Au. For layered materi- 
als having the badly conducting axis along the wire thickness 
(dashed red), only a slight improvement is gained, and the de- 
pendence on the anisotropy is negligible. If the badly conduct- 
ing axis is along the width of the wire (dashed-dotted blue), 
the improvement is more pronounced, but still saturates at rel- 
atively low anisotropy. For quasi-lD materials (dotted green - 
quasi-lD with both low conductivity terms of the same order; 
solid black - the more extreme case having a^^ <^ Oyy <C a^^x, 
where we assumed o^^ fa Oyy/r), the suppression is much more 
significant. For high anisotropy the scaling is linear with the 
anisotropy. Points: Examples of specific materials, not normal- 
ized to Au (see Table 



The suppression is somewhat more pronounced in the ex- 
treme case Gzz ^ CTyy *C <Txx- 

The reduction of the decoherence and heating rates is 
illustrated in Fig. [2] for the four cases discussed above. 
We plot the spin decoherence rate ([5]) for a super- 
position state of the hyperfine levels \F — 2,mp — 2) 
and \F = 2,mp = 1) in the ground state of ^^Rb. The 
anisotropic conductor is chosen such that the largest con- 
ductivity value axx coincides with the one for Au, an 
isotropic conductor taken as reference. The calculated 
rates for some specific materials are also given. It can be 
seen that quasi-lD materials are much more appealing to 
suppress heating and decoherence, although about one or- 
der of magnitude can already be gained with layered mate- 
rials, even at a relatively small anisotropy. In addition, for 
most anisotropic materials, even the best direction con- 
ducts worse than Au; therefore the rates presented in this 
graph even for layered materials are smaller than for Au. 

To conclude, heating and decoherence may be sup- 
pressed by several orders of magnitude even at room tem- 
perature, by using electrically anisotropic materials for 
current carrying structures on atom chips. 




Fig. 3. Geometry of the current carrying wire. Planar conduc- 
tivity fiuctuations can arise from bulk conductivity changes or 
thickness variations. A single Fourier component {kx,ky) is il- 
lustrated (color scheme represents conductivity changes) with 
the resultant variation of the current fiow direction. The cur- 
rent is tilted toward low resistivity regions and perpendicu- 
lar to the high resistivity ones. The conductivity fiuctuation 
is characterized by an angle 6, with tan(6') = kx/ky. Mea- 
surement and calculation for the isotropic case found that the 
largest amplitude of transverse currents is formed along wave 
fronts oriented a,t 9 — 45° . An intuitive explanation of this 
angle has been given in [34], and a formal model is developed 
in [36]. 

3 Time independent spatial corrugations of 
the magnetic potential 

Anisotropic materials affect not only the magnetic noise, 
leading to the time-dependent perturbations discussed so 
far, but also the static corrugations of the trapping po- 
tential that fragment an atom cloud. In this section, we 
discuss the current flow in imperfect wires, described by 
a spatially modulated conductivity and generalize the for- 
malism developed in Refs. [51U5S] to the anisotropic case. 



3.1 System definition 

We now consider a thin metal wire with conductivity fluc- 
tuations, as sketched in Fig. [3] The coordinate system is 
again defined such that the wire length L is along the x 
direction, its width W along y, and its thickness H along 
z. The anisotropic crystal used for the atom chip is again 
assumed to be 'aligned' with the wire axis, as in Eq. pUj) . 
This is of course not the most general scenario, however 
it is sufficient for our purposes here. 

3.2 Current corrugations 

In the isotropic case [34l[36] we start from Ohm's law 

E = pj, (22) 
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where p = a^^ = po + Sp is the scalar resistivity comprised 
of a homogeneous component po and a fluctuation Sp. This 
is plugged into Maxwell's equation, 



(23) 



We use a spectral expansion of the fluctuation 5p{x) into 
plane waves 

Ma;) = ^c"=-<5p(fc), (24) 

k 

where k = {kx, ky,kz) is a three-dimensional wave vector 
in which {ky,kz) take the discrete values 27r {ni / W, n / H) 
with integers m and n. Combining with the continuity 
equation V • J = we obtain to first order in Sp 



SJ{k) 



Spjk) 

PO 



Jo, 



(25) 



per each fc component. Here Jqx is the unperturbed cur- 
rent density applied along the wire. In the following, we fo- 
cus on thin (flat) wires and neglect the k^ component |36] . 
For a fluctuation 5p(k) with a given in-plane wave vector 
k = k {cos 6 , sin 6) , Eq. ([25|) implies a current fluctuation 
transverse to the applied current that can be characterized 
by the ratio 



a 



^(k) 



= isin(20) 



Sp(k) 



Po 



(26) 



For small perturbations, a describes the angular ampli- 
tude of the current deviation with respect to the x-axis. 
The main feature is that this deviation is maximal when 
the wave vector k is oriented at 9 = ±45°, as illustrated 
in Fig. [21 

The current density pattern is mapped onto the mag- 
netic field measured by the atoms at a distance d as 



1 



Sp{k) 



/3'"°(k,d) = e-'^''-sin(26') 

2 Po 



(27) 



exhibiting the same angular dependence. This is observed 
in experiments by imaging the atoms' density profile. The 
local column density profile is a direct measurement of 
the variations in the bottom of the trapping potential. 
The exponential term in /3 acts as a filter such that short 
wavelengths compared to the atom-surface distance are 
suppressed in the magnetic field fluctuations spectrum. 

In the case of electrically anisotropic materials, the re- 
sistivity becomes a tensor pi^i = pa.ki ^^Pki that is inverse 
to Eq. pU]) . Maxwell's equation combined with Ohm's law 
then takes the form 



[V X [pJ]]^ eijkdj [pkiJi) = 0, 



(28) 



where e^-fc is the Levi-Civita tensor and repeated indices 
are summed over. 

Rearranging and writing J = Jqx for the unperturbed 
applied current density, we obtain 



jkPo.kidjJi — eikj {djSpkx) Jo- 



(29) 



Again focusing on the 'aligned' case as in Eq. pO|) . and 
combining with the continuity equation diJi = 0, we solve 
the set of equations to get the components of the current 
fluctuations, 



'(fc) 



kx q 

fc • q 



Spxjk) 

PO,x 



Jo, 



(30) 



where q = {kx/ Po.x,kyl po,y,kz/ po,z) and 5px{k)/po,x is 
the relative fluctuation of the resistivity along the x- 
direction. This reduces to the isotropic result (P5)l when 

PO.x ~ PO,y ~ PO,z = Po- 

Focusing again on wave vectors k = k (cos 9, sin 9) and 
a thin wire, we get 



126* 



Spx 



r cos^ 



' Po.,: 



(31) 



where the anisotropy ratio is r = po,y/po,x- Here it is 
evident that the orientation of the current wave fronts 
is different from the isotropic case. While in the isotropic 
case, r = 1, we recover the sin(26') dependence of Eq. ((26l) . 
the angular dependence will asymptote to tan(0)/r in the 
highly anisotropic limit, r ^ 1. 



3.3 Analysis 

Let us introduce the angle ^max as that corresponding to 
the wave vector fc producing the largest current deviation. 
A simple calculation shows that this preferred angle sat- 
isfies tan6'inax = r^^^- In the highly anisotropic case (r 3> 
1), the current corrugation thus has wave fronts essentially 
perpendicular to the applied current (Fig.[4l lower panels). 
In addition, the peak corrugation amplitude, taken at the 
preferred angle and normalized to the isotropic case, scales 
like a''"'''°(fc,6lmax)/a''°(fc,45°) oc r'^/^. However, in an 
experiment, one rather observes the angle-averaged power 
spectrum for which we find a scaling 



ik,0)\ 



1/2 



Jd9\a'^° {k,9)\ 



oc r 



-3/4 



(32) 



in the high anisotropy limit r ^ 1. Highly anisotropic 
materials for trapping wires hence have the advantages of 
controlling the current scattering pattern in the wire and 
of suppressing fragmentation quite strongly. 

In Fig. 2] we show simulations of the fluctuations of the 
magnetic potential /3aniso {x^y^z = d) = e~'"*a™''°(fc), in 
arbitrary units. The plots are made in false color and 
consider different anisotropy ratios. We assumed a white 
spectrum for dp{k). The simulated wires are of width 
W = 200pm, thickness H = 2pm, and the observation 
point (atom-surface distance) is at d = 3.5pm. The scan is 
across the central lOOxGSOpm^, safely away from the wire 
edge. The change in orientation is clearly seen when the 
anisotropy ratio is varied. Note also the numerical values 
of the field corrugation amplitude (color bar), from which 
the overall suppression of fragmentation can be inferred. 
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Fig. 4. Comparison of magnetic corrugation patterns for dif- 
ferent anisotropy ratios r = axx/o'yy. The electron flow is 
mapped onto the magnetic field, shown here in false color 
(arb. units). Red colors indicate large magnetic field rota- 
tion /3™=° {x,y,z = d) = e-'=''a"-"'=°(fe) (Eq. (EIJ), and hence 
large transverse currents. As the anisotropy is larger, the 
value of the preferred angle Smax is shifted, following the law 
tanSmax = r^^^, and the overall signal intensity is suppressed, 
scaling on average as r~^^*. For the opposite case, r <^ 1 (badly 
conducting axis along the wire), the situation is reversed, and 
the fragmentation is enhanced. 



4 Materials and experimental design 
considerations 

In this section we discuss issues pertaining to the exper- 
imental verification of the theory, as well as possible ap- 
plications to quantum technology devices. We start with 
a review of electrically anisotropic materials which may 
be considered as potential candidates for atom chips, and 
then consider the effect of applying current to the wire, 
and hence heating it, on thermal noise properties. We con- 
clude by discussing the modifications imposed by the need 
to add longitudinal confinement to the trapping guide, in 
the form of a z-shaped wire. 



4.1 Electrically anisotropic materials 

Electrically anisotropic materials have been studied 
mostly in the context of characterizing their electrical 
transport properties (e.g. as a function of temperature 
or fabrication methods and parameters), or their mag- 
netic properties. A large portion of these materials are also 
high- Tc superconductors, hence they are also interesting 
in that context. Table 1 shows the relevant properties of 
some anisotropic materials. 

Although some hexagonal metals exhibit electrical 
anisotropy, it is usually smaller than r — A. Higher 
anisotropy ratios can be found in several layered com- 
pounds, as well as in some ladder-spin compounds, also 



shown in the tabic. A potential problem using these com- 
pounds is a strong dependence of electrical and magnetic 
properties on the exact relative composition, which could 
result in spatial inhomogeneity of the material properties 
within the wire. Paramagnetic materials are also consid- 
ered less preferable for magnetic traps. It should be noted 
however, that atom traps with permanent magnets have 
also been studied |67] . Some materials may have problems 
being in an ultra-high vacuum (UHV) environment, either 
due to out-gassing, or to a degradation of their electrical 
properties as the oxygen content in the ambient atmo- 
sphere is reduced (e.g. YBCO [68]). 

We note here especially different types of graphite, 
which seem to be attractive candidates for atom chips. 
Natural single crystal graphite |61] , exhibits an anisotropy 
of r « 100 and is non-magnetic. Highly ordered pyro- 
graphite (HOPG) is a form of synthetic single crystal 
graphite, which exhibits a very high (layered) anisotropy 
of r > 3500. This comes without significantly hindering 
of the good conductance along the wire, which is compa- 
rable for example to that of Ti, only one order of magni- 
tude lower than Au or Cu. However, the combination of 
graphite being a layered material with its softness creates 
a challenge in fabricating wire-sized structures out of it. 
Cutting the material perpendicular to the layers' plane 
may result in breaking the material into its layers (typical 
layer thickness on the order of ~ 30 nm for HOPG [64]). 
This has importance given that the anisotropy is perpen- 
dicular to the layers' plane. 

In Table 1 we also list as examples two materials with 
quasi-lD conductance. The Pcrovskite-typc transition- 
metal oxides SrNb03.4i, LaTi03.4i present relatively high 
conductivity in one direction, only one or two orders of 
magnitude below regular metals and much better than 
typical semiconductors. These are highly quasi-lD mate- 
rials, as in fact CTzz ^ o'yy <C (Jxx- These materials are 
commonly fabricated by the zone-melting method which 
enables an extremely high precision in the material com- 
position. This precision is required as materials with close 
structural proximity but different oxygen content have 
quite different electrical properties. In this method, high 
purity compound powders are combined at a high temper- 
ature through a scries of chemical reactions [69] ■ Special 
care is required to maintain exact weights of the various 
components and the desired oxygen content. The com- 
posite compounds arc then pressed and sintered into two 
rods. These are subjected together to a melting process, 
limited to a small region, in which one rod acts as a feed 
and the other as a seed. The melting zone is slowly trans- 
lated, and following the solidification of the melted parts, 
high-quality crystals can be obtained. For extensive re- 
views on these materials, including fabrication aspects, see 
[69ll70] . The common size of the crystals formed by this 
method is large compared with the requirements for atom 
chip micro-structures. However the electrical properties 
are maintained even upon size reduction by means of ap- 
plying pressure [7T][72] and various cutting methods. Still, 
in order to obtain the requirements for micro-structures 
more delicate methods, such as ion-beam etching, may be 
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Material Class 


Material 


Pa 


Pc 
Pa 


Ph 


Pa 
Ph 


Magnetic properties 


Source 




So 


27.6 


0.37 






Pfl rfl tn a OTi pt i c 


l47l 


11C|J iiieLaih 




/ .0 


91 


7 07 


9 9 


ryiama^iie Lie 






Te 


7nnn 


O.U4: 






l^iailldtiliC Lie 


1^01 




LaSb2 


65 


16.55 






Paramagnet ic 


|49l 


Layered compounds 


Sr^Ru907 
NaCo2 04 


1 on A 
ioU-40U 

3.8 

86 


ou - oUU 
23.5 
42.11 






Paramagnetic 


[DUIIDlllDZIIDo] 

[54] 
[55] 


Ladder-spin compounds 


Cai4 CU24 O41 

Sr3CaiiCu24 04i 


9-10'' 
4500 


0.01 
0.1 


10"^ 


10"* 


Paramagnetic 


[li 

[57] 




Bi2212 


180 


1.5 • 10* 






Paramagnetic 


EH] 


Cuprates 


YBCO 


90-360 


25-60 






Paramgnetic 


[59] 




La2-xSrxCu04 


45-400 


100 - 360 






Paramgnetic 


[60] 


Graphites 


Natural single crystal 


25.8 


105 






Nonmagnetic 


[61] 


HOPG 


18 


3750 






Nonmagnetic 




Perovskites 


LaTi08.4i 


362 


8.75 • 10* 


850 


lO""" 


Paramagnetic 


[65] 


SrNb03.4i 


226 


3- 10^ 


50 


2 • 10"^ 


Diamagnetic 


[66] 



Table 1. Electrically anisotropic materials as potential candidates for atom chip wires. The indices a, b, and c denote the 
crystal axes. We assume the good conductivity to be along the wire, and the worst to be perpendicular to the wire axis, in the 
plane parallel to the substrate (see definitions in Sect. 12.1]) . 



needed. These materials could be contacted electrically by 
standard means. To conclude, there seem to be no inher- 
ent hmitations in obtaining structures adequate for atom 
chips. 

4.2 Influence of wire heating on noise measurements 

An important practical issue that may compromise the re- 
duction of thermal noise is the heating of wires when cur- 
rent is applied. In principle, to investigate thermal noise it- 
self, no current is needed in the probed structure; however, 
considering an actual application with a magnetic micro- 
trap, there will be current in the wires, and hence heating 
of the surface. As can be seen from the prefactor of the 
magnetic field power spectrum (|20p . this heating would 
increase the noise, approximately linearly in the surface 
temperature. Bearing in mind that most anisotropic ma- 
terials have a lower conductivity than metals, this implies 
that running a current through an anisotropic material 
wire would hinder the improved lifetimes and coherence 
times. In Fig. [5] we plot as an example the expected life- 
time as a function of current density Jq along the wire, 
taking into account wire heating according to the model 
developed by Groth et al. [73j. As the relevant compo- 
nent to the heating process is the applied current, there 
is no difference between an anisotropic material and an 
isotropic one with the same conductivity along the wire, 
up to a factor 2 due to the anisotropy (Sect. 12. 3p . The 
longer lifetime relative to the Au case is mostly due to the 
lower conductivity relative to Au (in the direction along 
the wire for anisotropic materials), up to current densi- 
ties of lO^A/m^. For higher current densities we see that 
materials with low conductivity heat up much more than 
metallic ones (shown in the inset), and the correspond- 
ing lifetime drops by two orders of magnitude for current 
densities of up to lO^^A/m^, which is, however, already 
a fairly high value for small metallic wires [73]. It should 



be noted that the prefactor including the surface temper- 
ature in Eq. ([20]l in the context of trap loss is the same 
for atoms' heating- and decoherence-rates, hence this plot 
is relevant also for these processes. 

The conclusion here is that although electrically 
anisotropic materials are promising candidates for reduc- 
ing the coupling of noise to the atoms, their advantages 
are limited to moderate current densities. This should not 
pose overly stringent restrictions on possible applications, 
as at small atom-surface distances the required currents 
for small and tight traps are not very high. 



4.3 Contribution from structures providing longitudinal 
confinement 

The most common loffe-Pritchard magnetic micro-trap 
configuration is that of a z-shaped wire combined with an 
external bias field [Tj . The central part of the wire creates 
a guide potential, while the 'legs' provide the longitudinal 
confinement as well as the non-zero field value at the trap 
minimum, as illustrated in Fig. [5] In realistic magnetic 
micro-traps, where longitudinal-confinement is required, 
considering only the potential generated by the central 
part of the z-wire does not suffice for the anisotropic case, 
as the contribution from the 'legs' of the z-wire has to be 
taken into account. This is expected to be more and more 
important as the miniaturization of the micro-structures 
advances toward high packing density needed for exam- 
ple for matter-wave quantum technology devices. To asses 
the influence of the 'legs', we first consider a relatively 
large wire (atom-surface distance d = 50^m, wire width 
W = 100/J,TO, thickness H = 1/j,to, central part length 
Lcentrai = Imm, and 'legs' length Liegs = 1.5mm), repre- 
senting a typical micro-trap. To demonstrate the scenario 
of highly dense traps, we also consider the geometry of a 
much smaller structure (d — 3fim, W = 6/iTO, H ~ l/im. 
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Jo [A/m^] 

Fig. 5. Drawback of electrically anisotropic materials: reduced 
magnetic trap lifetime as a current-carrying wire heats up. The 
trap lifetime due to thermal noise induced spin flips is plotted 
as a function of current density Jq. For Jo — > 0, the level of 
pure thermal noise (at T = 300 K) is reached. The lifetime 
of metal wires (Au, solid red) is not significantly affected up 
to Jo ~ 10^^ A/rn^ . For high resistivity materials relative to 
Au the lifetime drops significantly as the current density is 
increased above lO^A/m^, due to wire heating. The lifetime of 
anisotropic materials (SrNb03.4i, solid black) is improved by a 
factor ~2, and their current density dependence is the same as 
that of isotropic materials of equal conductivity along the wire 
(dashed blue), (Sect. Inset: Temperature rise (above 300 
K) of the wire as the external current density is increased. Wire 
parameters are W = 10/i?n, H — 2.15/im, and trap distance 
d = This calculation follows the model developed by 

Groth et al. [73]. The substrate was assumed to be made of 
Si, with a 20 nm thick insulation layer of Si02. Steady state 
of the temperature increase was assumed to be reached after 
t = 30s. The resistivity was assumed to be independent of 
temperature in the AT range shown here, up to 150° above 
room-temperature [231166] . 



L central — 20 fim, and Liegs = 30^m). For both geome- 
tries we compare the spin decoherence rate as a function 
of material of the central part of the wire, normalized to 
that of a Au structure. We also take into account a current 
density Jq applied to the wire, limiting the heating of the 
wire to AT « 100°C {Sect.E^. 

The results are summarized in Table 2. Obviously it 
would be preferable to have the entire wire made of the 
anisotropic materials. However, for that the orientation of 
the crystalline axis in the area forming the 'legs' should be 
perpendicular to the central part of the z-wire, increasing 
fabrication complexity. However, having the 'legs' made 
of a normal metal, with isotropic conductivity, simplifies 
the fabrication, lowers power consumption and wire heat- 
ing. Indeed we find that even if the metal 'legs' to the 
anisotropic material reach up to the central part of the 
trapping wire, the reduction of the the decoherence rate 
is still significant. 




Anisotropic material 



Fig. 6. Schematic illustration of the studied z-shaped wires. 
The atoms are trapped at a distance d from the central part of 
the wire of length Lcentrai , made of an anisotropic material as 
summarized in Table 2. This part of the wire provides the radial 
trapping potential. The longitudinal confinement is provided 
by the z-wire 'legs' of length Liegs, made of Au. The width W 
and thickness H axe the same for the whole wire. 



Wire Material 


Large Wire 


Small Wire 


Central part 


Legs 


^spin decoherence 


'Yspin decoherence 


Au 


Au 


1 


1 


SrNbOi 


Au 


3.68-10"* 


9.9-10"^ 


SrNbOii 


Au 


4.39-10~'* 


1-10"^ 


HOPG 


Au 


9.7-10"'' 


1.8-10"^ 



Table 2. The affect of longitudinal confinement on the reduc- 
tion of the decoherence-rate when using electrically anisotropic 
materials. Two z-shaped wire geometries, as illustrated in Fig. 
[6] are considered: a large wire {d — 50/im, W = lOO/xm, 
H = l/im, Lcentrai ~ 1mm, and Li^gs ~ 1.5mm), and a small 
wire (d = S/xm, W = 6fj,m, H = lfj,m, Lcentrai = 20/im, 
Licgs = 30/im). The good conductivity axis (the a-axis) is 
along the central part of the wire. The spin decoherence rate, 
normalized to that of a structure entirely made of Au, is calcu- 
lated for different combinations of wire material compositions. 
The externally applied current density was set to 3.3-lO^A/m^ 
and 8.5 - lO^A/m^ for the large and small wires, respectively, 
such that the heating of the wire would not exceed 100° C (see 
section [42} . The difference between the two calculations for 
SrNb03.4i is whether the worst conductivity axis of this quasi- 
ID material is along the wire width (SrNbOi) or its thick- 
ness (SrNbOn). Even for small structures, with the Au 'legs' 
very close to the trapping area over the central part of the z- 
wire, the suppression of the decoherence and heating rates is 
still significant, further showing the appeal of using electrically 
anisotropic materials for realistic tight traps. 



5 Conclusions 

We have shown that utilizing electrically anisotropic 
materials on an atom chip will significantly suppress 
harmful mechanisms due to atom-surface interaction. 
Heating- and decoherence-rates due to thermal noise from 
the surface are expected to be strongly reduced even 
at room-temperature, scaling linearly with the electri- 
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cal anisotropy for certain types of anisotropic materi- 
als. This is maintained also in the presence of an ex- 
ternal current applied to the wire (causing wire heat- 
ing), up to current densities on the order of lO^A/rn^. 
We have shown that structures perpendicular to the 
wire axis, needed for longitudinal confinement, do not 
significantly hinder the improvements in heating and 
decoherence-rates, even when made of normal metal. The 
expected trap lifetime is not substantially modified by the 
use of such materials, although a significant improvement 
can be achieved by cooling the surface to low tempera- 
tures. Potential corrugations due to time independent cur- 
rent scattering patterns in the imperfect current carrying 
wires is expected to be suppressed with the anisotropy 
r as r~^/^ on average. The preferred orientation of scat- 
tered current wave fronts can also be controlled over a 
wide range by varying the electrical anisotropy, becom- 
ing perpendicular to the applied current in the extreme 
anisotropic case. Specific materials have been suggested as 
potential candidates for implementing the results obtained 
in this work. This could lead to improved functionality of 
atom chip based applications such as clocks, sensors and 
quantum information communication and processing. 
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